The Development of the SARS-CoV-2 Epidemic in Different Regions of Siberia in the 2020-2022 Period.

The comparison of the development of the SARS-CoV-2 epidemic in several neighboring regions can help researchers to assess the risks and develop more effective strategies and approaches in the field of preventive medicine. We analyzed the infection and mortality statistics for the 2020-2022 period in ten individual regions of the Siberian Federal District of Russia. We also sequenced complete genomes, which allowed us to analyze the genetic diversity of SARS-CoV-2 circulated in each of the ten regions and to build a phylogenetic dendrogram for the virus variants. The ParSeq v.1.0 software was developed to automate and speed up the processing and analysis of viral genomes. At the beginning of the pandemic, in the first two waves, the B.1.1 variant (20B) dominated in all regions of the Siberian Federal District. The third and fourth waves were caused by the Delta variant. Mortality during this period was at a maximum; the incidence was quite high, but the number of deposited genomes with GISAID during this period was extremely low. The maximum incidence was at the beginning of 2022, which corresponds to the arrival of the Omicron variant in the region. The BA.5.2 variant became the dominant one. In addition, by using NextClade, we identified three recombinants in the most densely populated areas.

Hydrodynamic description of protein folding: the decrease of the probability fluxes as an indicator of transition states in two-state folders.

Using hydrodynamic description of protein folding, the process of the first-passage folding of ubiquitin has been studied. Since a large number of folding trajectories were required to obtain converged folding flows, a coarse-grained representation of the protein in the form of a C-bead Go-model was employed, and discrete molecular dynamics was used to perform simulations. It has been found that the free energy surface has a maximum width in the transition state region, so that the densities of folding flows (probability fluxes) decrease to minimum when the system passes through the transition state. There are indications that the increasing number of different protein conformations in the transition state region compared with those in the neighboring regions of semi-compact and native-like states is responsible for the present phenomena. It has also been shown that if the free energy is projected onto a single reaction coordinate, the low populations of the transition states can be compensated by the increasing number of states, which can lead to a considerable decrease or even disappearance of the free energy barrier in the transition state.

Hydrodynamic description of protein folding.

A hydrodynamic description of protein folding is proposed and illustrated with a lattice protein model, which has a free energy surface (FES) typical of proteins with two-state folding kinetics. The flows from the unfolded to the native state are concentrated in a limited region of the FES. The rest is occupied by a flow "vortex", which does not lead to the native state. In contrast with intermediates that are associated with local minima, the vortex is not visible on the FES. The hydrodynamic interpretation thus provides new insights into the mechanism of protein folding and can be a useful complement to standard analyses.

A lattice protein with an amyloidogenic latent state: stability and folding kinetics.

We have designed a model lattice protein that has two stable folded states, the lower free energy native state and a latent state of somewhat higher energy. The two states have a sizable part of their structures in common (two "alpha-helices") and differ in the content of "alpha-helices" and "beta-strands" in the rest of their structures; i.e. for the native state, this part is alpha-helical, and for the latent state it is composed of beta-strands. Thus, the lattice protein free energy surface mimics that of amyloidogenic proteins that form well organized fibrils under appropriate conditions. A Go-like potential was used and the folding process was simulated with a Monte Carlo method. To gain insight into the equilibrium free energy surface and the folding kinetics, we have combined standard approaches (reduced free energy surfaces, contact maps, time-dependent populations of the characteristic states, and folding time distributions) with a new approach. The latter is based on a principal coordinate analysis of the entire set of contacts, which makes possible the introduction of unbiased reaction coordinates and the construction of a kinetic network for the folding process. The system is found to have four characteristic basins, namely a semicompact globule, an on-pathway intermediate (the bifurcation basin), and the native and latent states. The bifurcation basin is shallow and consists of the structure common to the native and latent states, with the rest disorganized. On the basis of the simulation results, a simple kinetic model describing the transitions between the characteristic states was developed, and the rate constants for the essential transitions were estimated. During the folding process the system dwells in the bifurcation basin for a relatively short time before it proceeds to the native or latent state. We suggest that such a bifurcation may occur generally for proteins in which native and latent states have a sizable part of their structures in common. Moreover, there is the possibility of introducing changes in the system (e.g., mutations), which guide the system toward the native or misfolded state.